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Abstract 

A search for pair-produced charged Higgs bosons is performed with the L3 detec- 
tor at LEP using data collected at centre-of-mass energies between 189 and 209 GeV, 
corresponding to an integrated luminosity of 629.4 pb~^. Decays into a charm and 
a strange quark or into a tau lepton and its neutrino are considered. No significant 
excess is observed and lower limits on the mass of the charged Higgs boson are 
derived at the 95% confidence level. They vary from 76.5 to 82.7 GeV, as a function 
of the H^ — »• TV branching ratio. 
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Introduction 



In the Standard Model of the electroweak interactions ^ the masses of bosons and fermions 
are explained by the Higgs mechanism 0. This implies the existence of one doublet of complex 
scalar fields which, in turn, leads to a single neutral scalar Higgs boson. To date, this Higgs 
boson has not been directly observed jHlll]- Some extensions to the Standard Model contain 
more than one Higgs doublet [Sj, and predict Higgs bosons which can be lighter than the 
Standard Model one and accessible at LEP. In particular, models with two complex Higgs 
doublets predict two charged Higgs bosons, H='=, which can be pair-produced in e"'"e~ collisions. 

The charged Higgs boson is expected to decay through H+ cs or H+ r+z/^^\ with a 
branching ratio which is a free parameter of the models. The process e^e~ — * H+H^ gives then 
rise to three different signatures: cscs, csr~z/^ and T^UrT^Ur- These experimental signatures 
have to be disentangled from the large background of the e^e~ — W^W~ process, characterised 
by similar final states. 

Data collected at centre-of-mass energies ^/s = 189 — 209 GeV are analysed here, supersed- 
ing previous results p. Data from i/i = 130 — 183 GeV jlj are included to obtain the final 
results. Results from other LEP experiments are given in Reference |H1 

The analyses do not depend of flavour tagging variables and are separately optimised for 
each of the three possible signatures. 

Data and Monte Carlo Samples 

The search for pair-produced charged Higgs bosons is performed using 629.4 pb~^ of data col- 
lected in the years from 1998 to 2000 with the L3 detector [^j at LEP, at several average 
centre-of-mass energies, detailed in Tabled 

The charged Higgs cross section is calculated using the HZHA Monte Carlo program jTIIj. As 
an example, at a/s = 206 GeV it varies from 0.28 pb for a Higgs mass, mH± , of 70 GeV to 0.17 pb 
for mH± = 80 GeV. To optimise selections and calculate efficiencies, samples of e"'"e^ H+H" 
events are generated with the PYTHIA Monte Carlo program for mH± between 50 and 
100 GeV, in steps of 5 GeV, and between 75 and 80 GeV, in steps of 1 GeV. About 1000 events 
for each final state are generated at each Higgs mass. For background studies, the following 
Monte Carlo generators are used: KK2f [12] for e"'"e~ — > qq(7), e"'"e~ — > yU"^/i~(7) and e"'"e~ — 
r+r-(7), BHWIDE HB] for e+e" ^ e+e", PYTHIA for e+e" ZZ and e+e" ^ Ze+e", 
YFSWW mi for e+e" W+W" and PHOJET (Tij and DIAG36 |16j for hadron and lepton 
production in two-photon interactions, respectively. The L3 detector response is simulated 
using the GEANT program ^1 which takes into account the effects of energy loss, multiple 
scattering and showering in the detector. Time-dependent detector inefficiencies, as monitored 
during the data taking period, are included in the simulations. 

Data Analysis 

The analyses for all three final states are updated since our previous publications at lower 
centre-of-mass energies jHlIZI: the searches in the H+H" cscs and csT'i^r channels are based 
on a mass dependent likelihood interpretation of data samples selected jTH] for the studies of 

^^The inclusion of the charge conjugate reactions is implied throughout this Letter. 
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W pair-production, while a discriminant variable is introduced for the search in the T^VrT Vt 
channel. These analyses are described below. 

Search in the H"""!!" — > cscs channel 

The search in the H+H^ cscs channel proceeds from a selection of high multiplicity events 
with balanced transverse and longitudinal momenta and with a visible energy which is a large 
fraction of ^/s. These criteria reject events from low-multiplicity processes like lepton pair- 
production, events from two-photon interactions and pair-production of W bosons where at 
least one boson decays into leptons. The events are forced into four jets by means of the 
DURHAM algorithm jTHj and a neural network JH] discriminates between events which are 
compatible with a four-jet topology and those from the large cross section e"^e~ qq(7) 
process in which four-jet events originate from hard gluon radiation. The neural network 
inputs are the event spherocity, the energies and widths of the most and least energetic jets, 
the difference between the energies of the second and third most energetic jets, the minimum 
multiplicity of calorimetric clusters and charged tracks for any jet, the value of the 1/34 parameter 
of the DURHAM algorithm and the compatibility with energy-momentum conservation in e^e~ 
collisions. After a cut on the output of the neural network, two constrained fits are performed. 
The first four-constraint fit enforces energy and momentum conservation, modifying the jet 
energies and directions. The second five-constraint (5C) fit imposes the additional constraint 
of the production of two equal mass particles. Among the three possible jet pairings, the 
one is retained which is most compatible with this equal mass hypothesis. Events with a low 
probability for the fit hypotheses are removed from the sample and a total of 5156 events are 
observed in data while 5112 are expected from Standard Model processes. The corresponding 
signal efficiencies are between 70% and 80%, for mH± = 60 — 95 GeV. 

Likelihood variables [20] are built to discriminate four-jet events compatible with charged 
Higgs production from the dominating background from W pair-production. A different like- 
lihood is prepared for each simulated Monte Carlo sample corresponding to a different Higgs 
boson mass. Seven variables are included in the likelihoods: 

• the minimum opening angle between paired jets; 

• the difference between the largest and smallest jet energies; 

• the difference between the di-jet masses; 

• the output of the neural network for the selection of four-jet events; 

• the absolute value of the cosine of the polar angle of the thrust vector; 

• the cosine of the polar angle at which the positive charged^^ boson is produced; 

• the value of the quantity 21n|M|, where M is the matrix element for the e"'"e~ 
W"'"W^ f our f ermions process from the EXCALIBUR [22] Monte Carlo program, 
calculated using the four-momenta of the reconstructed jets. 

Figures la, lb and Ic show the distributions of the last three variables while Figure Id 
presents the distribution of the hkelihood variable for mH± = 70 GeV. A cut at 0.7 on this 

Charge assignment is based on jet-charge techniques pi] . 
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variable, which maximizes the signal sensitivity, is applied as a final selection criterion, for 
all mass hypotheses. The numbers of observed and expected events are given in Table El and 
the selection efficiencies in Table 01 The main contributions to the background come from 
hadronic W-pair decays (70%) and from the e^e~ — > qq(7) process (26%). Figure 2 shows the 
5C mass of the pair-produced bosons before and after the cut on the final likelihoods. Peaks 
from pair-production of W as well as Z bosons are visible. 

Search in the H+H" — > csT~i>r channel 

The search in the H+H^ —>■ cst^Ut channel selects events with high multiplicity, two hadronic 
jets and a tau candidate. Tau candidates can be identified either as electrons or muons with 
momentum incompatible with that expected for leptons originating from direct semileptonic 
decay of W pairs, or with narrow, low multiplicity jets with at least one charged track, singled 
out from the hadronic background with a neural network ^Hj. The tau energy is reconstructed 
by imposing four-momentum conservation and enforcing the hypothesis of the production of two 
equal mass particles. The events must have a transverse missing momentum of at least 20 GeV 
and the absolute value of the cosine of the polar angle of the missing momentum is required to 
be less than 0.9. Finally, the di-jet invariant mass is required to be less than 100 GeV and the 
mass recoiling against the di-jet system less than 130 GeV, thus selecting 1026 events in data 
while 979 are expected from Standard Model processes, mainly from W pair-production where 
one of the W bosons decays into leptons and the other into hadrons. The signal efficiency is 
about 50%. 

To discriminate the signal from the background, mass dependent likelihoods [201 are built 
which contain eight variables: 

• the di-jet acoplanarity; 

• the angle of the tau fiight direction with respect to that of its parent boson in the rest 
frame of the latter; 

• the di-jet mass; 

• the quantity 2 In \M\ calculated using the four-momenta of the reconstructed jets and tau 
as well as the missing momentum and energy; 

• the transverse momentum of the event, normalised to ^/s^, 

• the polar angle of the hadronic system, multiphed by the charge of the reconstructed tau; 

• the sum of the angles between the tau candidate and the nearest jet and between the 
missing momentum and the nearest jet; 

• the energy of the tau candidate, calculated in the rest frame of its parent boson and scaled 
by y/s. 

The distributions of the last three variables are shown in Figures 3a, 3b and 3c. Figure 3d 
presents an example of the distributions of the likelihood variable for mH± = 70 GeV for data, 
background and signal Monte Carlo. A cut at 0.6 is applied for all likelihoods. This cut 
corresponds to the largest sensitivity to a charged Higgs signal. Table El gives the numbers of 
observed and expected events, while the selection efficiencies are given in Table El Over 95% 
of the background is due to W pair-production. Figure 4 shows the reconstructed mass of the 
pair-produced bosons before and after the cut on the final likelihoods. 
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Search in the H+H — > t^Ut'T channel 

The signature for the leptonic decay channel is a pair of tau leptons. These are identified either 
via their decay into electrons or muons, or as narrow jets. 

The selection criteria are similar to those used at lower jHlIZ]. Low multiplicity events 
with large missing energy and momentum are retained. To reduce lepton-pair background, an 
upper cut is placed on the value of the event collinearity angle, ^, defined as the maximum 
angle between any pair of tracks. The distribution of this variable is shown in Figure 5a. The 
contribution from cosmic muons is reduced by making use of information from the time-of-fiight 
system. Figure presents the distribution of the scaled visible energy, E^is/^/s, for events on 
which all other selection criteria are applied. 

The analysis is modified with respect to those previously published jHllZI in that the nor- 
malised transverse missing momentum of the event, Pt/Ems, whose distribution is shown in 
Figure Efc, is used as a linear discriminant variable on which no cut is applied. 

The efficiency of the H+H^ t^Ut-t^Ut selection for several Higgs masses is listed in Ta- 
ble 01 The numbers of observed and expected events are presented in Table El The background 
is mainly formed by W-pair production (60%), two-photon interactions (26%) and lepton pair- 
production (9%). 

Results 

The number of selected events in each decay channel is consistent with the number of events 
expected from Standard Model processes. A technique based on a log-likelihood ratio 4J is used 
to calculate a confidence level (CL) that the observed events are consistent with background 
expectations. For the cscs and csr~z/,- channels, the reconstructed mass distributions, shown 
in Figures Eb and Eb, are used in the calculation, whereas for the T~^i'T-T~Ur channel, the 
distribution of the normalised transverse missing momentum, shown in Figure Et, is used. 

The systematic uncertainties on the background level and the signal efficiencies are included 
in the confidence level calculation. These are due to finite Monte Carlo statistics and to the 
uncertainty on the background normalisation. The former uncertainty is 5% for the background 
and 2% for the signal Monte Carlo samples. The uncertainty on the background normalisation 
is 3% for the H+H^ cscs channel and 2% for the csr~i/^ and r^UrT^Ur channels. The 
systematic uncertainty on the signal efficiency due to the selection procedure is estimated by 
varying the selection criteria and is found to be less than 1%. These systematic uncertainties 
decrease the mH± sensitivity of the combined analysis by about 200 MeV. 

Figure IHl compares the resulting background confidence level, 1 — CL^, for the data to 
the expectation in the absence of a signal, for three values of the tv branching ratio: 

Br(H^ Tu) = 0, 0.5 and 1. The 68.3% and 95.4% probability bands expected in the absence 
of a signal are also displayed and denoted as la and 2a, respectively. A slight excess of 
data appears around mH± = 69 GeV for Br(H^ — > tu) = 0, as previously observed [H]. It is 
compatible with a 2.5cr upward fluctuation in the background. The excess is also compatible 
with a 2.9a downward fluctuation of the signal^-*. As observed in Figures 6b and 6c, no excess 
is present in the csT~Ur and t^v^-t^i)^ channels around mH± = 69 GeV. Therefore, the cscs 
excess is interpreted as a statistical fluctuation in the background and lower limits at the 95% 
CL on mH± are derived ^ as a function of Br(H^ — >■ rv). Data at ^/s = 130 — 183 GeV |7j are 
included to obtain the limits. Figure [7| shows the excluded mH± regions for each of the final 

■^^As an example, for Br(H^ tv) — 0.1, these figures are l.Scr and 2.7cr, respectively. 
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states and their combination, as a function of Br(H^ tu). Table 0] gives the observed and 
the median expected lower limits for several values of the branching ratio. 

In conclusion, refined analyses and larger centre-of-mass energies improve the sensitivity 
of the search for charged Higgs bosons produced in e"'"e~ collisions as compared to previous 
results |6i||7]. No significant excess is observed in data and a lower limit at 95% CL on the 
charged Higgs boson mass is obtained as 

mH± > 76.5 GeV, 

independent of its branching ratio. 
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(GeV) 


188.6 


191.6 


195.5 


199.5 


201.7 


204.9 


206.4 208.0 


Luminosity (pb~^) 


176.8 


29.8 


84.2 


83.3 


37.2 


79.0 


130.8 8.3 



Table 1: Average centre-of-mass energies and corresponding integrated luminosities. 





Channel 

CSCS CSr~Z/^ T~^UT-T~h'r 


Data 

Background 


2296 442 141 
2228 464 141 


Signal 


100 76 50 



Table 2: Number of observed data events and background expectations in the three analysis 
channels. The uncertainty on the background expectations is estimated to be 5%. The numbers 
of expected signal events for mH± = 70 GeV and Br(H''' — > ti/) — 0, 0.5 and 1 are also given 
for the CSCS, csr^P,- and t^Vt-t'Dt channels, respectively. 



Channel 



cscs 
csr~z/^ 



mH± 



Selection efficiency (%) 
60 GeV 70 GeV 80 GeV 90 GeV 95 GcV 



62 


62 


50 


58 


64 


38 


51 


43 


43 


39 


26 


30 


33 


34 


36 



Table 3: Selection efficiencies for various charged Higgs masses. The efficiencies are largely 
independent of the centre-of-mass energy. The uncertainty on each efficiency is estimated to 
be 2%. 



Br(H± ^ ru) 


Lower hmits (GeV) at 95% CL 
observed expected 


0.0 


76.7 


77.5 


0.26 


76.5 


75.6 


0.5 


76.6 


76.5 


1.0 


82.7 


84.6 



Table 4: Observed and expected lower limits at 95% CL for different values of the H='= — > tv 
branching ratio. The minimum observed limit is at Br(H^ tv) — 0.26. 
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2 In IMI Likelihood 



Figure 1: Distributions for the H+H~ — > cscs channel of: a) the absolute value of the cosine 
of the polar angle of the thrust axis, b) the cosine of the polar angle of the positively charged 
boson, c) the logarithm of the squared matrix element for the e"'"e~ W^W~ process and 
d) the selection likelihood for mH± = 70 GeV. The points represent the data and the open 
histogram the expected background. The hatched histogram indicates the expected distribution 
for a signal with mH± = 70 GeV and Br(H='= tv) = 0, multiphed by a factor of 10. The arrow- 
in d) shows the position of the cut. 
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Figure 2: Reconstructed mass spectra in the H+H~ — > cscs channel, for data and expected 
background, for events a) before, and b) after, the cut on the hkehhoods. The points represent 
the data and the open histogram the expected background. The expected distribution for 
mH± = 70 GeV and Br(H^ — > tu) = is shown as the hatched histogram. 
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Figure 3: Distribution for the H+H" — > csr~z/^ channel of: a) the cosine of the polar angle 
of the hadron system multiplied by the charge of the tau candidate, b) the sum of the angles 
between the tau candidate and the nearest jet and between the missing momentum and the 
nearest jet, c) the scaled energy of the tau candidate in the rest frame of the parent boson and 
d) the selection hkelihood for mH± = 70 GeV. The points represent the data and the open 
histogram the expected background. The hatched histogram indicates the expected distribution 
for mH± = 70 GeV and Br(H='= tu) = 0.5, multiplied by a factor of 5. The arrow in d) shows 
the position of the cut. 
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Figure 4: Reconstructed mass spectra in the H+H~ csr^P^ channel, for data and expected 
background, for events a) before, and b) after, the cut on the hkehhoods. The points represent 
the data and the open histogram the expected background. The expected distribution for 
mH± = 70 GeV and Br(H^ — > tu) = 0.5 is shown as the hatched histogram. 
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Figure 5: Distribution for the H+H^ t^Ut-t'Uj. channel of: a) the event colhnearity angle, ^, 
b) the scaled visible energy and c) the normalised transverse missing momentum of the event. In 
a) and b) all other selection criteria are applied and the arrows indicate the cut on the displayed 
variable. The points represent the data and the open histogram the expected background. 
The hatched histograms indicate the expected signal distributions for mH± = 70 GeV and 
Br(H± ^ Tu) - 1. 
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Figure 6: The background confidence level, 1 — CL^,, as a function of mH± for the data (solid 
line) and for the expectation in the absence of a signal (dashed line), for three values of the 
TV branching ratio. The shaded areas represent the symmetric la and 2a probability 
bands expected in the absence of a signal. 
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Figure 7: Excluded regions for the charged Higgs boson in the plane of the tv branching 

fraction versus mass, for the analyses of each final state and their combination. The dashed line 
indicates the median expected limit in the absence of a signal. Regions below mH± = 50 GeV 
are excluded by data collected at the Z resonance [SBj and at y/s = 130 — 183 GeV |2j. 
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